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Abstract 

r-| ■ Weakly interacting massive particles (WIMPs) remain a prime candidate for 

P-1 the cosmological dark matter (DM), even in the absence of current collider 

JIL[ signals that would unambiguously point to new physics below the TeV scale. 

The self-annihilation of these particles in astronomical targets may leave 
observable imprints in cosmic rays of various kinds. In this review, we focus 
on gamma rays which we argue to play a pronounced role among the various 
possible messengers. We discuss the most promising spectral and spatial 
signatures to look for, give an update on the current state of gamma-ray 
searches for DM and an outlook concerning future prospects. We also assess 
in some detail the implications of a potential signal identification for particle 
DM models as well as for our understanding of structure formation. Special 
emphasis is put on the possible evidence for a 130 GeV line-like signal that 
was recently identified in the data of the Fermi gamma-ray space telescope. 



1. Introduction 

Evidence for a sizable non-baryonic and cold dark matter (DM) compo- 
nent in the universe derives from an impressive range of unrelated cosmolog- 
ical observations [l| , covering distance scales from tens of kpc to several Gpc 
and leaving very little room for alternative explanations. On cosmological 
scales, DM contributes a fraction of VL-^ = 0.229 ± 0.015 to the total energy 
density of the universe Weakly Interacting Massive Particles (WIMPs) 
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provide a theoretically particularly appealing class of candidates for the so 
far obscure nature of DM j^, with the lightest supersymmetric neutralino 
often taken as a useful template for such a WIMP. It is often argued that 
the thermal production of WIMPs in the early universe generically leads to 
a relic density that coincides with the observed order of magnitude of Q^, 
though this rests on the assumption of a standard cosmological expansion 
history and there exist well-motivated particle physics scenarios that predict 
alternative production mechanisms for WIMP DM While the LHC non- 
observation of new particles below the TeV scale (apart, possibly, from the 
Higgs boson) has already prompted doubts whether the WIMP DM scenario 
is still our best bet [s], it must be stressed that electroweak low-energy ob- 
servables {g — 2 in particular) do favor new physics contributions not too far 
above 100 GeV 0] . While this tension starts to considerably disfavor very 
constrained models of, e.g., supersymmetry p], it may simply be an indica- 
tion that the new physics sector which the WIMP belongs to appears at a 
much smaller mass scale than any new colored sector. 

Attempts to identify WIMP DM can be classified into collider searches 
for missing transverse energy, direct searches for the recoil of WIMPs off the 
nuclei of terrestrial detectors and indirect methods that aim at spotting the 
products of WIMP self-annihilation. Among possible messengers for such 
indirect searches, gamma rays play a pronounced role as they propagate es- 
sentially unperturbed through the galaxy and therefore directly point to their 
sources, leading to distinctive spatial signatures; an even more important as- 
pect, as we will see, is the appearance of pronounced spectral signatures. 
This prime role of gamma rays provides our motivation for an updated and 
dedicated review on these messengers, which we hope will prove useful and 
complementary to existing general reviews on indirect DM searches In- 
deed, the recent indication for a DM signature in gamma-ray observations of 
the Galactic center (GC) j^, |lo[ makes such a review extremely timely, and 
we therefore dedicate a considerable part of it to discuss in great detail both 
the status of the potential signal and its implications. 

Gamma rays can either be observed directly from space or, via the showers 
of secondary particles they trigger in the atmosphere, indirectly with ground- 
based experiments. The former option necessarily implies rather small effec- 
tive areas and an upper bound on the photon energy that can reliably be 
resolved, but allows for a large field of view and the observation of gamma 
rays at comparably small energies. Particularly promising instruments for 
the latter option are imaging Air Cherenkov Telescopes (lACTs) that detect 
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Table 1: Rough comparison of basic telescope characteristics relevant for indirect DM 
searches with gamma rays, for a selection of typical space- and ground-based experiments 
that are currently operating, shortly upcoming or planned for the future. The quoted 
sensitivity is for point sources at the 5cr level, after lyr (50 hrs) of space- (ground-) 
based observations and assuming typical backgrounds. Where applicable, numbers refer 
to photon energies at or above E ~ 100 GeV (1 TeV). More details in Refs. |16| (Fermi- 
LAT), [i3 (AMS-02), [H (GAMMA-400), ^ (MAGIC), ^2^] (HESS-II) and [21| (CTA). 



the Cherenkov light emitted by the shower particles and use efficient image 
reconstruction algorithms to determine the characteristics of the primary 
photon. These instruments have a limited field of view and a lower energy 
threshold set by the need of discriminating photons from the background of 
primary muons and hadronic cosmic rays; their extremely large effective area 
and rather small field of view make them ideal for pointed observations. In 
Table [H we pick typical examples for space- and ground-based experiments 
that are currently operating or planned for the future and compare some 
basic telescope characteristics that are particularly relevant for DM searches. 
Experiments that fall into the same broad categories but are not listed ex- 



as 



plicitly in the Table include for example AGIL E and VERITAS [12!^, _ . 
well as the future CALET [lij and DAMPE [IJlIS- We stress that the 



numbers in Tab. [T]are intended to provide a convenient order-of-magnitude 
comparison of instrumental characteristics; they should not be used as the 
basis of detailed sensitivity estimates (see, however, the stated references). 

The expected DM-induced gamma-ray fiux from a direction ip, averaged 
over the opening angle ^ip of the detector, is given by 

where the integration is performed along the line of sight (l.o.s.), (cri;)ann is the 
average velocity-weighted annihilation cross section, m-^ the mass of the DM 
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particle (for which we assume X = x) j Px density, Bf the branching 

ratio into channel / and A^;^ the number of photons per annihilation. An often 
quoted reference value for (o"f)ann is the so-called 'thermal cross section' of 
(av) ~ 3 ■ 10~^^cm^s~^, which is the annihilation rate expected for thermally 
produced WIMPs in the most simple case (i.e. s-wave annihilation without 
resonances or co-annihilations j22|). The right part (in parentheses) of the 
Eq. dl]) contains all the particle physics input and, for the typically very 
small DM velocities, is usually sufficiently independent of v{r) that it can be 
pulled outside the integrals (note, however, that this is not true for strongly 
velocity-dependent cross-sections like in the case of Sommerfeld enhancement 



23l |24| . resonances or p-wave annihilation). It contains the full spectral 
information that we will discuss in some detail in Section The remaining 
part, sometimes referred to as the astrophysical factor, contains in that case 
the full information about the spatial distribution of the signal and will be 
discussed in Section [31 

We continue by reviewing in Section H] gamma-ray limits on DM anni- 
hilation as well as the current status of claimed DM signals. The poten- 
tially enormous implications of a signal identification for our understanding 
of both structure formation and the underlying particle model are then out- 
lined in Section [5l with a focus on the intriguing 130 GeV feature in the 
direction of the GC. We discuss future prospects for the detection of DM 
with gamma rays in Section [H] and conclude in Section [71 For most of this 
review, we will assume that DM consists of WIMPs; many aspects, however, 
can be applied — or generalized in a straight-forward way — to other cases as 
well, most notably decaying DM jisi (for which one simply has to replace 
|(crf)p^ — 7- m^Tp^ in Eq. ([T|), where F is the decay rate). Where apphcable, 
we will comment on this on the way. 

2. Spectral signatures 

At tree level, DM particles annihilate into pairs of quarks, leptons, Higgs 
and weak gauge bosons. The hadronization and further decay of these pri- 
mary annihilation products leads to the appearance of secondary photons, 
mainly through 7r° — >■ 77, and the resulting gamma-ray spectrum dNj^/dE^ 
can be obtained from event generators like Pythia Codes like DarkSUSY 



27j provide user-friendly numerical interpolations of these spectra, based on 
a large number of Pythia runs, but there also exist several analytic param- 
eterizations in the literature jl, [isj . Secondary photons show a featureless 
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spectrum with a rather soft cutoff at the kinematical hmit = and are 
universal in the sense that dNj^/dx (with x = E^/m-^ takes a very similar 
form for almost all channels / and only very weakly depends on m^. A con- 
vincing claim of DM detection based exclusively on this signal, which would 
show up as a broad bump-like excess over the often rather poorly understood 
astrophysical background, appears generically rather challenging. 

For this reason, it is often much better warranted to focus on the pro- 
nounced spectral features that are additionally expected in many DM models— 
not only because they greatly help to discriminate signals from backgrounds, 
and hence effectively increase the sensitivity of gamma-ray telescopes to DM 



signals [29|, but also because a detection may reveal a lot about the under- 
lying model for the particle nature of the DM. Before discussing the vari- 
ous types of spectral features in more detail, let us briefly mention further 
contributions to the total photon yield that do not give rise to such nice 
spectral features but may still visibly change the spectrum (in particular 
at E^ <^ ^x)- ■'■^ models with large branching fractions into e^e~ pairs, 
e.g., inverse Compton scattering of those on starlight and the cosmic mi- 



crowave background leads to additional gamma rays [30[. Another source of 



additional low-energy photons are electroweak [31| and strong [32| radiative 
corrections, with an additional gauge boson in the final state; these contri- 
butions can actually be quite sizable for DM masses much larger than the 
gauge boson mass, in particular if the tree-level annihilation into a pair of 
SM particles is suppressed. 

2.1. Lines 

The direct annihilation of DM pairs into jX — where X = Z,H or 
some new neutral state — leads to monochromatic gamma rays with E^ = 

[l — / , providing a striking signature which is essentially im- 
possible to mimic by astrophysical contributions |33[. Unfortunately, these 
processes are loop-suppressed with 0{al^) and thus usually subdominant, 
i.e. not actually visible against the continuous (both astrophysical and DM 
induced) background when taking into account realistic detector resolutions; 



however, examples of particularly strong line signals exist [2J, |3J|. A space- 
based detector with resolution AE/E = 0.1 (0.01) could, e.g., start to dis- 
criminate between 77 and 7Z lines for DM masses of roughly < 150 GeV 
[m^ < 400 GeV) if at least one of the lines has a statistical significance of 



> 5(7 [35|. This would, in principle, open the fascinating possibility of doing 



'DM spectroscopy' (see also Section EJ. 
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2.2. Internal hremsstrahlung (IB) 

Whenever DM annihilates into charged particles, additional final state 
photons appear at (9(aem) that generically dominate the spectrum at high 
energies. One may distinguish between final state radiation (FSR) and virtual 



internal hremsstrahlung (VIB) in a gauge-invariant way j36|, where the latter 
can very loosely be associated to photons radiated from charged virtual parti- 
cles. FSR is dominated by coUinear photons, thus most pronounced for light 
final state particles, m/ ^ m^, andproduces a model-independent spectrum 
with a sharp cut-off at = [37|, |38|; a typical example for a spectrum 



dominated by these contributions is Kaluza-Klein DM [39i]. VIB, on the 
other hand, dominates if the tree-level annihilation rate is supp ressed (like 
e.g. the annihilation of Majorana particles into light fermions |40|) and/or the 
final state consists of bosons and the t-channel particle is almost degenerate 



with [4l|. It generates pronounced bump-like features at < which 



closely resemble a slightly distorted line for energy resolutions AE/E > 0.1. 



The exact form of VIB spectra, however, is rather model-dependent [36 — 
which in principle would allow an efficient discrimination between DM models 
for large enough statistics (see e.g. (42|). 

2.3. Cascade decays 

Another possibility to produce pronounced spectral features is DM anni- 
hilating into intermediate neutral states, XX ~^ which then decay directly 



((/)-)■ 77 j43|) or via FSR (e.g. (p -)■ £+£"7 j44|) into photons. While the 
latter situation results in a spectrum that resembles the standard FSR case 
(with a slightly less sharp cut-off and a potentially considerably reduced rate 
in the degenerate case, ~ m^), the former process induces a box-shaped 



spectrum with a width of AE = yrn'^ — m|; for small mass differences, it is 
thus indistinguishable from a line. 

In Fig. [H we compare the various different spectra discussed above; in 
order not to overload this figure, however, we do not include the FSR- 



dominated spectrum from lepton final states (see e.g. Ref. j39|). Secondary 
photon spectra from all possible quark or weak gauge boson final states are all 
contained in the rather thin gray band (we adopted = 100 GeV, though 
the result is quite insensitive to this value). For the VIB spectrum, we as- 
sumed Majorana DM annihilation into light fermions via a scalar t-channel 
particle ('sfermion') almost degenerate in mass with x, like encountered in 
supersymmetry 36(|, and for the box we chose m<^ = 0.95m^. 



6 



Bringmann & Weniger (2012) 



AE/E = 0.15 



10 



AE/E = 0.02 




5 



qg-,ZZ,WW- 



0.01 



0.02 



0.05 



0.10 



0.20 



0.50 



1.00 



2.00 



X = E / m 



'X 



Figure 1: Various gamma-ray spectra expected from DM annihilation, all normalized to 
N{x > 0.1) = 1. Spectra from secondary particles (gray band) are hardly distinguishable. 
Pronounced peaks near the kinematical endpoint can have different origins, but detectors 
with very good energy resolutions AE/E may be needed to discriminate amongst them in 
the (typical) situation of limited statistics. See text for more details about these spectra. 

3. Spatial Signatures 

The peculiar morphology of annihilation signals, tracing directly the DM 
density, offers another convenient handle for discriminating signals from back- 
grounds. The most relevant targets are the GC, dwarf spheroidal galaxies 
and galaxy clusters with respective half light radii of roughly 61/2 ^ 10°, 
61/2 ~ 0.1° and 61/2 ^ 0.1°. Further important targets are DM clumps or the 
angular power spectrum of the isotropic gamma-ray background (IGRB), all 
of which we will discuss in this section. 

3.1. Halo Profiles and the Galactic Center 

The arguably brightest source of gamma rays from DM annihilation is the 
center of our Galaxy. Within a few degrees around the GC, WlMPs would 
induce a gamma-ray flux of about (9(10""^) ph cm~^ s~^ at the Earth (at 
> 1 GeV, assuming a thermal annihilation rate into bb, = 100 GeV and 
standard halo profiles), very well in reach of current instruments. However, 
the line-of-sight to the GC traverses the galactic disc, which harbours nu- 
merous high-energetic processes (tt^ production in cosmic-ray interactions. 
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Bremsstrahlung and inverse Compton emission, bright point sources); the 
corresponding gamma-ray fluxes of (9(10~^) ph cm~^ s~^ thus outshines DM 
signals in many cases by orders of magnitude. Furthermore, the uncertainties 
in the signal and background morphologies make the identification of a DM 
signal from the inner Galaxy a challenging task. 

A useful general parametrization of DM halos, which encompasses a large 
number of commonly used profiles, reads 



r 


-7 


'l + (rQ/r,)" 






1 + {r/r-sY 



(2) 



where r is the distance from the halo center, Tq ^ 8.5 kpc the position of 
the Sun and ~ 0.4 GeV cm~^ the local DM density j45| (see Ref. 46| for 
a discussion of systematic uncertainties of this quantity). The parameters a, 
P and 7 determine the halo shape, and the concentration. The commonly 
used Navarro- Frenk- White (NEW) profile 47|], for example, is obtained for 
(a,/3,7) = (1,3, 1) (with ~ 20 kpc in case of the Milky Way); the cored 
isothermal profile follows when setting (a,/3,7) = (2,2,0) and ~ 3.5 kpc 
(note, however, that matching observational constraints in principle results 
in a rather large range of allowed values for jisj ). Typically, kinematic 
observations of line-of-sight velocities do not sufficiently constrain the DM 
profile, so one has to rely on (extrapolated) results from numerical A^-b ody 
simulations of dissipationless structure formation in a ACDM cosmology 49|, 



50 ; see also Ref. 51 for a review. Most recent results tend to somewhat 



favor a spherical Einasto profile [52 



Einasto^^^ 



P@ exp 



^OlE 



aE 



^CtE 



(3) 



over the somewhat steeper NEW profile. Both and the halo shape param- 
eter ue depend on the total halo mass; in case_of the Milky Way, numerical 
simulations yield ue — 0.17 and ~ 20 kpc [49]. Note that for r < O.Olr^, 
there are actually numerical indications for profiles with logarithmic slopes 
of about 7 ~ 1.2, i.e. steeper than both NEW and Einasto [49|, |53 . 



In the Milky Way, baryonic matter dominates the gravitational potential 



for roughly r 



< 



05 



which can have a great impact on the DM distribution 



with respect to the expectations from DM-only simulations mentioned above. 
In particular, the cooling and infall of baryons could — by a mechanism known 
as adiabatic contraction — lead to a steepening of the inner DM profile 54 . 
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However, such a scenario becomes much less hkely if, as sometimes found in 
simulations 55|, feedback from star formation and supernovae dominates over 
cooling and infall processes. In fact, the presence of baryons could have the 
opposite effect of producing cores rather than cusps, see e.g. the discussion 
in Ref. 56 1, even for a system with the size of the Milky Way js^] (though 
such a core might also form only in the very center of a contracted profile 
[isl). The adiabatic growth of the central supermassive black hole (SMBH) 
could also generate a central spike of DM within the inner ~ 10 pc if 
the SMBH seed starts out close to the GC |t60i] (SMBH mergers, on the other 
hand, would rather destroy initial cusps in the profile |6l|)- Microlensing 
and stellar rotation curve observations can only exclude the most extreme 
scenarios, yielding upper limits of about 7 < 1.5 for the logarithmic slope of 
the DM profile near the GC j62|. The resulting difference in the expected 
DM annihilation flux from the inner ~ 0.1° around the GC, when comparing 
the most extreme cases of a cored profile and a profile as steep as this upper 
limit, amounts to around five orders of magnitude 28|]. In any case, the 
annihilation signal from the GC would most likely appear as an extended 
source with a peculiar angular profile 63|. Due to the large astrophysical 



foreground in the very center (e.g. the bright HESS source J1745-290 64|). 
the optimal region of interest (ROI) for signal searches extends out to a few 
degrees and could also lie slightly away from the GC j§, 29, 65 



A topic that has received only little attention is the possibility that the 
point of highest DM density could be displaced from the GC. The latter lies 
at the dynamical center of our Galaxy, coinciding with the position of Sgr 
A*. Baryons are directly affected by star formation and supernovae, and 
they shock during galaxy mergers whereas DM does not. It is not unlikely 
that at least during some periods in the merger and formation history of our 
Galaxy significant displacements existed 59|], like e.g. observed for the SMBH 



in M87 If such a displacement remained until today, it could provide 
a spectacular window into the formation history of our Galaxy, but would 
also introduce another unknown into the search for a DM signal. Actually, 
state-of-the-art simulations of late-time spiral galaxies with a significant bar, 
like the Milky Way, show first hints for such a displacement 58, 67l . 



3.2. Substructure Enhancement 

For annihilating DM, the gamma-ray emissivity is proportional to the DM 
density squared. Unresolved substructures in the DM distribution, confirmed 
to exist by all cold DM A^-body simulations, can hence potentially have a 
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huge impact on the signal strength [68| - simply because the astrophysical 



factor in Eq. (JT]) effectively averages over and one always has (p^) > (p^)^ 
for inhomogeneous distributions. This effect is typically quantified in terms 
of the so-called boost factor B, which can be defined as the ratio of the actual 
line-of-sight integral to the one obtained assuming a smooth (e.g. Einasto) 
component only. 

Roughly, every decade in subhalo masses down to the cutoff scale con- 
tributes the same to the total gamma-ray flux [69| (with small subhalos 



possibly being slightly more important [51[), though the details depend cru- 
cially on the adopted subhalo mass distribution and concentration, as well as 
survival probabilities of the smallest clumps (all of which have to be extrapo- 
lated over many orders of magnitude from the results of A^-body simulations). 
The lower cutoff in the subhalo mass distribution is set by the kinetic de- 
coupling of WIMPs in the early Universe [TQI] and strongly depends on the 
DM particle properties | 7ll . 
scale of dwarf galaxies [73| 
for standard MSSM neutralinos 172 



while in principle it could be as large as the 
it falls into the range of roughly lO^^^-lO^^M© 



Not too much is known about the precise distribution of substructures, 
but what one can learn from A^-body simulations is that due to merger and 



tidal stripping in the Milky Way halo (see e.g. |50l. l74| and references therein), 
fewer substructures are expected in the inner galaxy than in the outer parts 
(though the surviving subhalos close to the halo center have larger concen- 
trations). This implies that the expected boost factor for GC observations 
is of order unity, while it may be as large as (9(1000) for galaxy clusters 
due to the enormous contained hierarchy of masses |75l-l77l| . Note also that 
the expected angular dependence of the signal can change in the presence of 
large boost factors: in the limit where unresolved substructures completely 
dominate the total signal, the flux essentially scales with (the line-of-sight 
integral over) rather than p^; for the Milky Way, this implies that the 



halo flux emissivity could be changed for r > 1 kpc [7^, 178|, [79 



3.3. Point- like Sources 

Many complications associated with the GC are avoided when looking at 
point-like targets outside the galactic disk. The corresponding signals are 
typically considerably fainter, which is however compensated by the greatly 
simplified and much smaller astrophysical background. 

The probably most promising source class are nearby dwarf spheroidal 
galaxies. These faint satellites of the Milky Way exhibit the largest known 



10 



mass-to-light ratios, up to ~ 1000 Mq/Lq, and do not show signs for gas 
or recent star formation. As such, they are not expected to be gamma- 
ray emitters 80|. A canonical set of less than ten dwarf spheroidals were 
subject of numerous recent studies [8l|, |82|. Since dwarf spheroidals are DM 



dominated^ stellar kinematics can be efficiently used to constrain the DM 
Despite the remaining uncertainties in the shape of the DM 



content 



83, 84 



profile, it turns out that the integrated signal fluxes are surprisingly robust 
(at the 10%-50% level) js^, 85| (though a substructure boost may enhance 
the signal by up to two orders of magnitude in the most optimistic case 



This is in sharp contrast to the large uncertainties related to the signal flux 
from the GC that we discussed above and makes dwarf spheroidal galaxies 
excellent targets to derive robust constraints on the DM annihilation cross- 
section. Furthermore, it allows a simultaneous analysis of multiple dwarf 
galaxies, which further increases the sensitivity for DM signals. 

Galaxy clusters are the most massive DM dominated virialized objects 
in the Universe and provide excellent targets to search for an annihilation 



signal [75|, l87| . Harboring an enormous hierarchy of substructures, they are 
the astronomical targets that are expected to maximize the boost factor 
In optimistic scenarios 77| , they could therefore outshine a signal from 



77 



local dwarf spheroidals by a factor of a few [85|, |88|, which makes them 



very attractive as targets for the potential detection of a DM signal [75 



77|. The rather large involved astrophysical uncertainties connected to both 
the subhalo distribution and cosmic-ray induced gamma rays, on the other 
hand, imply that robust limits derived from cluster observations are usually 
not competitive. Similar to dwarfs, the sensitivity to cluster signals can be 
enhanced by a combined analysis SSNOO . 

Given that they come with no intrinsic astrophysical backgrounds, clumps 
(subhalos) of DM in the galactic halo that are not massive enough to trigger 



star formation are further important targets for indirect DM searches |91 



If discovered as unidentifled sources with no counterparts at other wave- 
lengths by surveying instruments like the Fermi-LvVT [92|] , detailed follow-up 
observations with lACTs could become vital to prove their DM nature by 
means of additional angular and spectral information. 

3.4- Extragalactic Diffuse Signal and Anisotropics 

Gamma rays from DM annihilation at cosmological distances, integrated 
over all redshifts and (sub) halo distributions, appear largely isotropic and 
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add up to the astrophysical IGRB together with contributions from galac- 



tic DM annihilation |93l-l96| . The astrophysical IGRB is presently not very 
well understood, but believed to stem from unresolved sources like blazars, 
star-forming galaxies and milli-second pulsars. A WIMP with thermal an- 
nihilation cross-section and (9(100 GeV) masses could contribute between 
~ 1% and ~ 100% to the IGRB, depending on the minimum mass and 
abundance of subhalos j97| (though the intrinsic uncertainties related to the 
substructure distribution may well be even larger ^8|). In case of annihila- 
tion into 77, the signal would be a gamma-ray line broadened by the redshift 
and provide a peculiar spectral signature to look for in the IGRB |93|, which 
can be efficiently constrained by observations |97 . 

An interesting approach towa rds signal identification is also to exploit its 
angular power spectrum 98l - ll0l| , which receives contributions from subhalos 
within our own Galaxy as well as from extragalactic (sub)halos. Let us 
mention here in particular that anisotropy measurenaents might turn out to 
be a feasible way to probe the minimal subhalo mass 100l | , which would open 



a completely complementary window into the particle nature of the DM [72 



4. Current status 
4.I. Limits 

The total number of photons above the detector threshold, typically dom- 
inated by secondary photons, is a very convenient and simple measure to con- 
strain possible exotic contributions to observed gamma-ray fluxes. Limits on 
the DM annihilation rate are therefore usually presented in the (av) vs. 
plane, with the assumption of WIMPs dominantly annihilating into bb being 
an often adopted standard that is useful for comparison. Such c onstraint s 



have been derived fr om the ob servation of galax y clusters |88l . l89l . Il02l . Il03 



or the IGRB [96|, l97|, Ill9l . 



external galaxies |28l. Il04 
lite dwarf galaxie s |82l. |85 



120 



105i. globular cluste rs llOdllOTi . Milky Way satel 



108Hm|, the GC |ll2l - lll6j and halo [117|, [118 



The currently best limits of this kind, for WIMP masses 5 GeV < < 
1 TeV, derive from observations of nearby dwarf galaxies by the Fermi satel- 
lite [ssQ for m^^ < 25 GeV, these limits are actually stronger than the 'ther- 



mal' rate of {o'v)^i ~ 3 ■ 10 



At ~ 700 GeV, they weaken to 



^ Depending on the assumed profile and small-scale cutoff ( as w ell as subhalo prop- 
erties), the recently presented constraints from galaxy clusters |l03| are nominally even 
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{av),-, < 4 ■ 10 



and for even higher WIMP masses, the currently 



strongest hmits a re pr esented by the HESS collaboration from observations 
of the GC region |ll5| : at ~ 1 TeV (10 TeV), thos e are abo ut a factor of 



10 (30) weaker than the thermal value (see also Ref. (109|, Il2ll | for indepen- 
dent studies of these limits). When comparing limits from different targets, 
however, one should always keep in mind that the underlying astrophysical 
uncertainties that enter as the line-of-sight integral in Eq. ([T]) may be quite 
different; in particular, as stressed in Section |3l predictions for integrated 
signal fluxes are much more robust for dwarf galaxies than for the GC. 
There have also been various sear ches for line signals: in M31 with 



HEGRA 122 , at the GC with EGRET 123 , and with Fermi-LAT GC 124 



126| as well as dwarf data |127j | and in galaxy clusters 128|. The currently 



strongest limits presented by t he LAT collaboration follow from Fermi obser- 



vations of the GC region 126| and extend from (cxu)^^ < 3 



10 



-2() S - 

cm s 



at 



= 10 GeV to (aw)^^ < 4 - 10 



-27 S - 

cm s 



at 



200 GeV (slightly stronger 
limits can be found in independent analyses [10|, |125| for masses from 1 to 



300 GeV). Preliminary results from HESS exclud e line s above 500 GeV down 
to cross-sections of {(Jv)^^ ~ 2 x 10 



-27 S 

cm s 



1291. So far, none of those 



limits gets close to the expectation for vanilla WIMP models; realistic models 
featuring particularly strong line signals, however, start to get constrained. 

The small-scale {i > 150) gamma-ray anisotropics observed by Fermi- 
LAT, indicating the presence of some unresolved sou rce p opulation, can only 
partially (if at all) be explained by DM annihilation 120l | ; on the other hand , 



this can already be used to constrain the distribution of DM subhalos [101 



While some unidentified gar ama- ray sources in principle qualify as candidates 
for annihilating DM clumps 130|, the presence of any unconventional sources 
in current gamma-ray data seems very unlikely after taking into account 



surveys at other wavelengths [92]. Finally, let us mention that the assumed 



non-observation of gamma-ray point sources from DM annihilation limits 



the allowed abundance of ultracompact mini halos [131| . which can be used 



to put extremely stringent constraints on the power of primordial density 



perturbations [132 



tighter. There are also claims that limits from GC [113| or globular cluster [l07l | observa- 
tions are actually stronger (though seemingly much less robust) than the dwarf limits. 
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4-2. Signals 

Historically, there have been a couple of claims of potential DM signals in 
gamma rays. Fo r the GC , e.g. , they correspond to DM mas ses in the ~MeV 



1331, -lOGeV [113|,11MI1351, -100 GeV [1361, ~TeV 



times with possible counterparts at radio frequencies |l40l ) or even multi-TeV 



39l. I137H1391 (some- 



range 



141 



Possible DM signals in gamma rays have al so b een claimed in 
the diffuse gamma-ray flux jl42 | or from galaxy clusters jl43 |. While it can 



be argued that there still remains a bit of controversy in some of these cases, 
evidence is certainly not compelling; instead, in the past, more refined anal- 
yses and ne w data tended to disfavor the DM hypotheses previously put 
forward 6J, Il03l Il44j4l47l (for a di scuss ion of recent DM signal claims not 



only in gamma rays, see also Ref. 1481). O'^^ reason for this is that the 



claimed signals typically rely on the presence of some broad excess in the 
differential gamma-ray flux that was mostly attributed to secondary photons 
which, as stressed before, makes the identification of a DM signal intrinsi- 
cally error-prone |§ In this context, it is also worth recalling that neither very 
small {m^ <C 100 GeV) nor very large (m^ ^ 1 TeV) DM masses are eas- 
ily accommodated in realistic WIMP frameworks that successfully address 
shortcomings of the standard model of par ticle physics. Furthermore, one 
shoul d ap preciate the fact that both CMB 15 ll and cosmic ray antiproton 
data 152l | provide very stringent constraints on the possibility of (9(10) GeV 
WIMP DM — even though such a possibility might be interesting from the 



point of view of direct DM searches |153 



^mma-ray signal at 
9, QOl, on the other 



The recently discovered hint for a monochromatic 
around 130 GeV in the Fermi data of the GC region 
hand, would correspond to a rather natural DM mass of 0(100) GeV and, 
even more importantly, for the first time provide evidence for a gamma- 
ray feature which is widely regarded as a smoking gun signature for DM 
[33|. Performing a spectral shape analysis in target regions close to the GC 
(selected in a data-driven approach by using photons at much lower energies), 
the signal was found to correspond to a DM mass of = 149 ± 4 GeV 
for an assumed VIB signal jsj and = 129.8 ± 2.4 GeV for a 77 line 



lOj . in each case with a local (global) significance of almost 5 a (more than 



^ The only exception is the 511 keV hne f rom e~^e~ annihilation seen by Integral |149| . 
Its observed non-spherical distribution |l46| . however, makes a DM interpretation highly 
unlikely — which in any case would be restricted to a very narrow mass range, TOp < Wy S , 
3 MeV, in order not to overproduce continuum photons from final state radiation |37l . Il50l | . 
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3 a). The deduced annihilation rate depends on the DM profile; for an 
Einasto profile, e.g., it is {av)i+i~^ = (5.2 ± 1.3 x 10^^'^ cm^ s"^ and 



e.g. 

(av)^^ = (1.27±0.32 toil) x 10 
Tab. |2]for a dis cussi on of 'jZ and 'jH final states). This excess was confirmed 
independently 154j by adopting a different statistical technique based on 
kernel smoothing; this analysis also demonstrated that the intrinsic signal 
width cannot be much larger than the energy resolution of Fermi LAT (~10% 
at 130 GeV) — leaving only a VIB signal, a gamma-ray line or a narrow box 
as po ssible explanation in terms of DM annihilation. Later, Su &: Finkbeiner 



-27 pj^3 g 1^ respectively (see Section I5l2] and 



155| adopted a refined spatial template analysis to demonstrate that the 
existence of a gamma-ray line emitting region of radius ~3° close to the GC 
is preferred over the no line hypothesis with a global significance of more 
than 5(7 (under the assumption of an Einasto profile). 

It is worth emphasizing that the above described signal is the only signif- 
icant line-like feature in the sky that we currently find distinguishable in the 
data from around 20 GeV to at least 300 GeV (we checked this explicitly by 
performing line searches along the galactic disc, as well as a subsamphng anal- 
ysis of anti-GC data — see also Refs. 0, [l^, 155| : conflicting claims 154j . 15d | 
may likely be explained as statistical fluctuations at the expected level). 
However, it is quite interesting to note that there might be weak evidence 
for line signals, with much smaller significance but at the same energy, in 



the direction of galaxy clusters [l57l |. Further weak evidence for such lines 
has also been foun d in some of the unassociated gamma-ray point sources 
obse rved by Fermi 



158 



159 



These indications are currently intensely debated 
if eventually confirmed with better statistics than currently available, 
they would significantly strengthen a DM interpretatio n. Let us also mention 
that there is no c orrelation between the Fermi bubbles 
signal at the GC 148l Il55| ; the significant overlap 



139. 160 and the line 



161 of the bubbles with 



the target regions adopted in Refs. j^l, [l^] is thus purely accidental and re- 
lated to the peculiar angular distribution of signal and background photons 



fsee also Ref. 162 



As already stressed, the intrinsic signal width is small: assuming a Gaus- 
sian instead of a monochromatic signal, we find an upper limit of 18% at 
95%CL (though a pair of lines might provide a marginally better fit, see 
Section [5l2|) . In Fig. |2l we dem onstr ate furthermore explicitly that a broken 



power-law, as suggested in Ref. |161| , does not provide a reasonable fit to the 
data unless one allows for an absurdly large spectral break A7 = 72—71 ^ 10; 
the best fit is obtained in the line-like limit A7 — )■ 00. As also indicated in 
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Figure 2: This plot shows the confidence contours obtained when fitting the 130 GeV 
signature with a broken power-law with spectral break from 71 to 72 at 130 GeV (plus the 
usual background power-law with free normalization and slope) . Best fits are obtained for 
A7 = 72 — 7i ^ 10, when the signal approaches a line-like shape; the gray area indicates 
parameters that are realistically accessible for astrophysical sources. 



the figure, the required spectral indices both above (72) and below (71) the 
break would in that case be well outside the range of values expected from 
standard astrophysical processes (the power- law background in our fit, on 
the other hand, has a spectral index consistent with ~ 2.6, i.e. the value 
expected from cosmic-ray proton collisions with the interstellar medium). |f| 
Let us now briefiy turn to what currently appears as the greatest chal- 
lenge to a DM interpretation of the observed signal (of course, this necessari ly 
refiects our own bias; for extended discussions see e.g. Refs. |9l.Eol. [l5i . lT63) l 
The first caveat is that the signal center seems to be displ aced, by ab out 1.5° 
or 200 pc, from the dynamical center of the galaxy [155) (see Ref. (Tel 

for 

the claim that this could merely be a statistical artifact); even if such a dis- 
placement might in principle be possible (see discussion in Section IXT]) . this 



•^To generate the plot, we redid the analysis from Ref. [13] in Reg4 (SOURCE events), 
replacing the monochromatic line by a broken power-law that changes its spectral index 
from 71 to 72 at 130 GeV. Note that a smooth change of 7 would make the fit quality 
even worse and similar conclusions hold for the more typical case of a power-law with an 
exponential cutoff (which roughly corresponds to 72 1 in our parameterization). Fits 
are performed in an energy window from 60 to 250 GeV. 
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certainly came as a surprise. An even stronger threat to the DM hypothesis 
might be the indication, so far at a weaker level of significance, for a line 
in part of the gamma rays from cosmic-ray induced air showers in Earth' s 



atmosphere (commonly referred to as Earth limb or Earth albedo) |155l . |163 



However, the problematic limb photons only appear at a very specific range, 
30° ^ 6 ^ 45°, of incidence angles (unlike the signal from the GC which 
shows up at all 6). Furthermore, the majority of events with these incidence 
angles come actually not from the limb but are of astrophysical origin and 
this larger sample does not show any evidence for a 130 GeV feature ^163i] . 
This confusing situation might well be an indication that the limb excess is 
merely a statistical fluctuation that soon will disa ppea r with more limb data. 

So far no compelling alternative instrumental 165| or astrophysical mech- 
anism has been proposed that could actually produce such a line at 130 GeV 



see Ref. |166l | for an interesting proof-of-principle with fine-tuned pulsar 



winds - which however cannot explain the extended morphology of the sig- 
nal). We stress that all analyses of the line signal so far rely on the publicly 
available Fermi data and information only, and that line searches operate by 
construction at the statistical and systematical limitations of the instrument. 
There has not yet been an official statement from the Fermi collaboration 
concerning the signal, in particular with respect to whether the energy re- 
construction of Pass 7 events is reliable at energies above 100 GeV in light 
of the recent findings]^ Eventually, such an independent confirmation of the 
130 GeV excess will of course be indispensable. 



5. What could we learn from a signal? 

Gamma rays may carry important and nontrivial information about the 
nature of the DM particles. Let us now demonstrate in more detail what kind 
of information could actually be extracted in case of a signal identification. 
For definiteness, we will take the tentative line signal as an example and 
assume in this Section that it can indeed be explained by DM. 



^ Note that the latest official compilation of Fermi line limits 126l | was finalized before 
the first indication for the line signal was announced; it relies on 24 rather than 43 
months of data and excludes a significant part of the ROI close to the GC. The tentative 
signal claim is thus not in tension with those limits fld^ . 
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Figure 3: Comparison of different flux profiles as function of the opening angle A9 of 
an hourglass-shaped ROI that is centered on the GC (see text for detailed definition). In 
green we show the zblcr uncertainty band of the hne flux measured inside this region: while 
compatible with a standard Einasto profile at 6* > 1°, it is incompatible with both a cored 
and a sufhciently contracted profile, as well as with a signal from DM decay. 

5.1. Dark Matter distribution 

A gamma-ray line would allow to study the distribution of DM in the 
GC with unprecedented accuracy, which could serve as important feedback 
for state-of-the-art numerical simulations of gravitational clustering. To il- 
lustrate this, we show in Fig. |3]the zblcr range of the line flux as function of 
the opening angle ^ of a ROI centered on the GC (green band)0 Obviously, 
for large galactocentric distances, the flux drops because the GC signal is 
washed out. Remarkably, the flux profile is perfectly consistent with the 
predictions from a standard Einasto profile (red solid line). The dash-dotted 
black line shows for comparison the prediction for a DM annihilation signal 
from a sufficiently contracted profile (chosen to be essentially equivalent, for 
the angular resolution of Fermi-LAT, to a point source at the GC); at angles 
0^5° the measured fiux starts to deviate from the predictions, indicating 
that the 130 GeV signal is not a point source but extended up to these an- 



^ We use here an ROI with hourglass-like shape, defined hy ip < min(3°, 9) plus ip <9 
and \b/£\ > 0.7 (V' is the angular distance from the GC). Otherwise, we use the same 
procedure as in Ref. [13] to obtain the line fiux, i.e. a power-law + line fit to SOURCE 
class events. 
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gles. On the other hand, the signal is too concentrated to be compatible with 
a cored profile (dotted black line, here for a core radius of 3.5 kpc). Note 
that a contribution from Milky Way subhalos could further boost the signal 
at angles 9 > 20° by a factor of a few 7^, which however is not observed 
and thus might be used to place constraints on subhalo models. DM decay 
(blue dashed line) would also lead to a flux that is too weak at the GC to be 
compatibl e wit h the observations (unless extreme assumptions on the profile 



are made [l67l |. which however are likely to be in conflict with microlensing 
and dynamical constraints |162|). 

As already mentioned, there is s ome evidence for a ~ 1.5° offset of the 
signal with respect to the GC jl55 |. While this observation was certainly 
unexpected, it might possibly be explained by the interplay between baryons 
forming a bar and DM [67 1. T he same bar would however also likely destroy 



a DM cusp in the center |168| : tidal disruption by the central supermassive 
black hole may be a further issue to worry ^about. More data, as well as 

are thus needed to settle in 



58 



more detailed simulations along the lines of 
how far the morphology of the 130 GeV feature is compatible with theoretical 
expectations. In fact, such an improved understanding could eventually allow 
to infer important details about the formation history of our Galaxy. 

Once the GC signal is established, an exciting future apphcation would 
be a precise all-sky survey to look for the same 130 GeV feature, aiming 
at a (partial) map of the Galactic and cosmological DM distribution. For 
sufficiently large substructure boosts, the 130 GeV feature could for example 
appear as a bump in the IGRB or in galaxy clusters. So far, no corresponding 
lines were found in th e IGRB 97|] , but there is a possible weak indication from 
galaxy clusters 157l |. If confirmed, it would necessarily imply a rather small 
value for the small-scale cutoff in the subhalo distribution in order to produce 
the required large boost factors of 0(10^); this, in turn, could be used to 
obtain highly complementary bounds on the underlying WIMP model. With 
more data, it may also be possible to identify a few individual DM subhalos 
(see e.g. the heavily disputed we ak indica tions for a 130 GeV signal from 
unidentified sources of the 2FGL flSsl. IT59|l The galactic distribution and 
total number of those subhalos would provide invaluable information on the 
DM distribution and allow to discriminate between subhalo models that are 
currently discussed. 
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5.2. Dark Matter models 

At the time of this writing, the hterature has aheady seen a consider- 
able amount of model-building efforts to explain the line signal in terms 
of annihilating DM. This range s from p henomenological and in some sense 



model-independen t approac hes |ll4 Il69 



tive field theories 170-172 



or analyses in the context of effec- 
to concrete model building. Proposed solutions 
that mostly fall into the latter category include an additio nal U jl) symmetry 



173| . DM as t he lig htest state of a new scalar multiplet 



sneutrino DM 
DM 



175 , axion-mediated DM annihilation 



1741 ]. right-handed 
1761 1 . two- component 



1771 ]. magnetic inelastic DM |178l ] or dipole-interacting DM [179|. Even 



neutralinos have been proposed as a possible cause of the sign al, a lbeit in 



non-minir aal versions of supersymmetry like no-scale J-'-SU(5) |180| or the 



NMSSM |l7ll . 1 1761 . Il8lj — each time, however, arguing for additional indi- 
cations in favor of the respective model in collider data. The possibility 
of d e caying D M being responsible for the signal has also been entertained 



1671 . Il82l Il83l ] — though the expected angular dependence of the signal in this 
case is hardly consistent with observations, see the previous subsection. 

One generic problem for any realistic model-building is that the annihila- 
tion cross section required to fit the data is considerably larger than typically 
expected for thermally produced DM, at least if the rehc density is set by the 
tree-level annihilation rate. On the particle physics side, possible ways to en- 
hance the annihilation rate in that case include the Sommerfeld enhancement 



23l |24| in the presence of new light bosonic messenger particles that mediate 



an attractive force between the initial state DM particles or, at the cost of 
some fine- tuning, the presence of a resonance (i.e. s-channel annihilation 
via a new neutral particle with m ~ 2m^, the same spin and CP properties 
as th e initial state); yet another mechanism might be cascade annihilation 
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On the astrophysical side, larger annihilation fluxes arise by adopting 
a larger local DM density or a profile that is steeper in the innermost part 
than our reference Einasto profile, Eq. ([3]); see, however. Fig. [3] for the rela- 
tively tight constraints on the latter option. Note that even if one relaxes the 
theoretically appealing assumption of thermal DM production, one needs to 
worry about large annihilation rates at tree-level because they would produce 
secondary photons potential ly in stark conflic t with co ntinuum g amm a-ray 



data 114, 116, 128 



1671 . Il85| ; also antiproton |l67l . Il86j and radio |186| data 



are quite efficient in constraining such large annihilation rates. This fact can 
be used to rule out e.g. Wino or Higgsino DM as an explanation for the line 
(see also below). Antiprotons could be constraining for future experiments 
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Figure 4: SUSY scan comparing the expected number of quasi-monochromatic photons 
(120 GeV < < 140 GeV) to the number of secondary photons. Green (red, blue) points 
correspond to models where thermal production leads to a reUc density in (smaller than, 
larger than) the observed range. Filled symbols indicate models where VIB contributes 
at least 3 times more photons than 77 and 7Z. Exclusion limits |ll6| and signal fit (at 
la and 2a) both assume an Einasto profile. Dashed lines show the effect of enhancing the 
annihilation flux by the same amount in both exclusion (|6|, |^| < 5°) and signal ROI; only 
models below the dashed line may thus in principle explain the line at 130 GeV. 



not only because of the tree-level annihilation rate, but also due to the asso- 
ciated DM annihilation into qq fi nal states if the 77 signal is dominated by 



colored particles in the loop |187 



In order to illustrate the ab ove p oint, we consider in Fig. Hlthe result of a 
large scan (for details, see Ref. 188[) over the parameter space of the cMSSM 
and a phenomenological MSSM-7, keeping only neutralino DM models where 
IB, 77 and 7Z photons for Ey G [120, 140] GeV dominate the secondary 
contribution by a factor of at least 5; for reference, we also show the case 
of pure Wino (W) and Higgsino ('-f^') DM. Filled symbols correspond to 
models where IB photons outnumber line photons by at least a factor of 
3. Assuming an Einasto profile as in Eq. (|3]), we also show the required 
signal stren gth to account for the line observation (shaded area) as well as 
limits 116| on the continuum fiux from the GC region (solid lines). Note 
that both limits and signal region are roughly proportional to (/ dQds p^) , 
albeit integrated over slightly different regions near the GC; assuming this 
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factor to be the same for both regions, dashed hnes indicate how hmits 
and signal region would change if adopting a profile that is different from 
Eq. ([3]). From this figure, we can draw at least three important conclusions: 
i) As already anticipated, the large annihilation rate required to explain the 
signal cannot easily be achieved for thermally produced neutralino DM. ii) 
Even when enhancing the annihilation rate such as to sufficiently increase 
the production of 77 or 7Z final states, it is very difficult to do so without 
violating the bounds from continuum gamma rays. In fact, the observed 
correlation between loop- and tree-level rates is generally expected from the 



optical theorem and should thus not only apply to neutralino DM |186| . iii) 



VIB, on the other hand, does not follow this pattern and can thus be argued 
to be a more natural explanation for such a strong line-like signal (still in 
need, however, of (9(10) enhancement factors for standarcl^ neutralino DM). 

The possibly only way to avoid the above considerations may be to 
strongly restrict any coupling of DM to charged standard model particles 



lighter than loop signals could then easily dominate over VIB 



signals without being in conflict with constraints arising from tree-level an- 
nihilations. In such a cas e, o ne would rather generically expect not only 
one but at least two lines |170| and the observed ratio of photon counts (or 
limits on th ose) can p rovide crucial information about the underlying par- 



ticle model 170, 191 



In Fig. [5] we therefore provide significance contours 
and upper limits for a second hne besides the observed 130 GeV feature; for 
convenience, we summarize these results in Tab. |2] in terms of limits on the 
annihilation cross section {av).yx under the assumption that the signal cor- 
responds to DM ami i hilat ion into jY (for X, F = 7, Z, h). Interestingly, as 
observed earlier 155 . 170l |. one can see a weak indication for a second line at 



^ Let us stress that all supersymmetric models of Fig. 2] assume unification of gauge cou- 
plings at the GUT scale, which prohibits a large Wino fraction of the lightest neutralino. 
Scanning a simple phenomenological MSSM-9 ^189.] , where this assumption is relaxed, we 
found the top right branch of non-thermally produced models (open red symbols in Fig. |4]) 
to extend all the way to the pure Wino case — albeit always above the (extended) dashed 
line. However, we could not find any models with larger IB rates than shown in Fig. |4] 
(possibly due to the restricted nature of the scan and MSSM version employed). 
Secondary photons from electroweak and strong corrections, see Section [2l were not in- 
cluded and would move the VIB dominated models in Fig. 2] somewhat upwards. While 
a dedicated future analysis is certainly warranted, let us stress that we still expect the 
continuum gamma-ray limits to be easily satisfied for thermally produced neutralinos 
(antiproton constraints from these channels [190.] are likely less stringent |9|]). 
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Figure 5: Significance contour (thick black lines in la steps) and upper limits (yellow line; 
95%CL) for a second line. Assuming that the 130 GeV feature is due to 77 (green sohd), 
7Z (blue dashed) or jH (red dotted), the vertical lines show the corresponding positions 
of the other two lines. A very weak hint for a ■yZ line at 114 GeV can be identified. 



7X 


[GeV] 


{av)^x [10-27cm=^s-i] 






X 


77 
7Z 
7i7 


129.8 ± 2.4+^4 
144.2 ± 2.21^2 
155.1 ±2.11^1 


1.27 ±0.32+°;^^ 
3.14±0.79l[J:^[! 
3.63 ± OMtlil 


1 

< 0.28 

< 0.17 


1 

< 0.79 


< 0.83 

< 1.08 
1 



Table 2: Upper limits at 95%CL (or best-fit value with ±lcr error) on the branching 
ratios into the secondary line, assuming that the primary line at E-y — 130 GeV is due to 
annihilation into 7X with X = Z 01 H. Note that {av) ^ z I {fyu) -y-y < 2.01 at 95%CL. 



114 GeV — which coincides surprisingly well with the energy expected for a 
7Z line if the 130 GeV feature can be attributed to DM annihilation into 77; 
for this case, we also state the best fit value for the ratio of cross sections. 

6. Future prospects 
6.1. Next decade 

The next ten years will bring a plethora of new results in indirect DM 
searches. It is right now that experiments start to probe vanilla WIMP DM 



23 



models and thus will either identify a signal or exclude many of the most com- 
mon scenarios. Ongoing experiments like Fermi-LAT, HESS-II, VERITAS 
and MAGIC will continue to take data, may identify new targets for DM 
searches, profit from a better understanding of astrophysical backgrounds 
and prepare the stage for upcoming instruments like CTA or GAMMA-400 
with considerably improved characteristics for DM searches. Indirect detec- 
tion with gamma rays will also profit from an interplay with upcoming results 
from neutrino searches with IceCube, anti-matter searches with AMS-02, re- 
sults from the LHC as well as from next-generation direct WIMP detectors. 
Furthermore, continuously improving results from N-body simulations that 
realistically take into account the various components of baryonic matter will 
sharpen our understanding of the signal morphology. 

Assuming a ten year lifetime of Fermi-LAT, the limits on the annihilation 
cross-section that were derived from observations of nearby dwarf galaxies 
with 2 years of data 
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would improve on purely statistical grounds by a 
factor of y/b to 5, depending on the annihilation channel and the DM mass 
(which determines whether the limits are derived in the signal-dominated 
high-energy regime or in the low-energy regime dominated by the diffuse 
gamma-ray background). Optical surveys like Pan-STARRS, the Dark En- 
ergy Survey or the Stromlo Missing Satellite Survey could increase the num- 
ber of known dwarf spheroidals by a factor of 3, which could additionally 
incre ase t he constraining power by a factor of y/S to 3 in the most optimistic 
case 
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Further significant improvements are expected from the upcoming 
Pass 8 version of the LAT event reconstruction software, which will lead to an 
enhanced effective area for high energy gamm a ray s due to a better hadron 



rejection and an improved energy resolution |193| . It is hence conceivable 
that Fermi-LAT dwarf limits will improve by a factor up to ten, which could 
allow to constrain WIMPs with thermal annihilation rate into bb up to DM 
masses of ~ 600 GeV. Similar improvements might be expec ted for limits 



from galaxy clusters |89l . Il02 



116. 185 and the halo 



DM searches in the GC 
on the other hand, will mostly profit from a refined understanding of 
astrophysical backgrounds. The best limits on annihilation into gamma-ray 
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lines 10 



or VIB features Q are right now based on almost four years of 
data; more data and improved event reconstruction will strengthen them by 
at least a factor of v^. 

By now, HESS observations of t he G C provide the strongest limits on DM 
annihilation with > 700 GeV 115| . The newly mounted 28m-diameter 
telescope HESS-II has an about 3 times lower energy threshold than HESS- 
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I, as well as ad dition al timing information which will improve the rejection 



of cosmic rays [194| . Both improvements will help to extend the current 
constraints to lower masses and values of the annihilation cross-section. If 
Fermi-LAT identifies a DM signal candidate at high enough energies, HESS 
could quickly confirm it thanks to its large effective area, and provide ad- 
ditional information about the variability and spatial extent of the source. 
Further results are anticipated from VERITAS, MAGIC and AMS-02 as well. 
A substantial improvement in the TeV regime will eventually come with 
CTA. For an NEW or Einasto profile, observations of the GC are expected 
to exclude cross-sections down to the thermal one at TeV DM masses — an 
improvement of m ore t han an order of magnitude with respect to the current 



HESS constraints [195| . Similar improvements with respect to HESS are also 
expected for dwarf galaxy observations, although they would still hardly be 
competitive with results from space-based instruments. 

Future space-based instruments like GAMMA-400 or CALET/DAMPE 
will — thanks to an extended imaging calorimeter and a large lever arm to 
the converter foils — have a much better energy and angular resolution than 
Fermi LAT. However, they will come with a somewhat smaller effective area. 
For that reason, DM limits from e.g. dwarf spheroidal observations would 
likely only be somewhat strengthened in the background limited regime at 
low energies. However, these instruments would be excellent machines for 
detailed follow up studies of DM signal candidates that might be identified 
in the Fermi-LAT data, in particular in case of pronounced spectral features. 

In Fig. |6] we provide a convenient summary plot of limits on the DM 
annihilation cross section into bb (in red and blue for space- and ground- 
based instruments, respectively) as well as in 77 (in green). The limits 
are collected as a representative selection of different instruments, and we 
concentrate on observations of dwarf spheroidal galaxies, the GC and the 
Galactic halo, since they provide right now the most stringent constraints. 
All limits are shown as function of the time of their publication, and were 
derived for common assumptions on the DM profile (NFW profiles in most 
cases); the different symbols correspond to limits for DM masses of 10, 100 
and 1000 GeV. Limits to the right of the dashed black line show expected 
limits for the next decade, as well as the fundamental reach to be discussed 
in Sec. 16.21 In the past eight years, most limits have improved by an order 
of magnitude; a similar improvement is expected during the upcoming ten 
years. Even without excessive boost factors, these limits start to reach deep 
into the parameter space of WIMP DM models. In particular observations 
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Fermi Halo (CPS) [iihFermi Halo (col.) [HI; Whipple Draco llOSll: Veritas Willman 



I 



mj : Veritas Segue I |l97j: CTA Segue I and GC HESS GC 



115|; EGRET 77 T23| 

Fermi 77 lyr [ijj; Fermi 77 2yr (VW) [l25|; Fermi 77 2yr (col.) |l26| : GAMMA-400 [351] 



of the GC with Cherenkov telescopes at high energies as well as observations 
of dwarf spheroidal galaxies with space-based instruments at lower energies 
have a great potential for deriving constraints or discovering a signal. 

Finally, the prospects for a further study of the 130 GeV feature in the 
Fermi-LAT data, if it persists, are extremely good. HESS-II has just seen 
its first light and given the good performance foreseen for the instrument 
in hybrid mode, it should allow a quick confirmation of the signal reported 

In 



m 



10 



if systematic uncertainties are sufficiently under control [35 
the case of the future space-based telescopes GAMMA-400 and DAMPE, the 
improved energy resolution provides an enormous potential not only for the 
detection of the 130 GeV feature, but also for the efficient discrimination of 
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a VIB feature from one or several lines 

A substantial substructure boost will generally be necessary even in the 
near fut ure t o measure the 130 GeV signal elsewhere in the sky, like in dwarf 
galaxies 127j , the EGBG 97| , galaxy clusters or Milky Way subhalos. O n the 



other hand, the recent claims for an identical signal in galaxy clusters |157 
or unidentified sources of the 2FGL 158l | are likely to soon be confirmed or 
refuted in light of the upcoming data and improved telescope performances 
mentioned above. 
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6.2. Fundamental reach 

Let us finally estimate what may be called the systematic or "fundamen- 
tal" reach of gamma-ray searches for DM signals. As an instr uctiv e exercise, 
and somewhat complementary to other works (see e.g. Ref. flSs) ). we will 



here assume infinite observational time (or effective area), and estimate the 
remaining systematic limitations. These limitations come from (1) a model- 
ing of astrophysical backgrounds, and (2) the instrument itself. 



Numerical codes like GALPROP |198l | do an excellent job in modeling 



the galactic diffuse gamma-ray emission. However, uncertainties related to 
e.g. the simplified propagation set-up, the interstellar radiation field or mod- 
els for the gas distribution inhibit a priori predictions, and even after per- 
formir ig fits to the data up to ~30% residuals at large and small scales re- 



mam 



1181 . |199| . At high latitudes, away from the galactic disk, background 
variations at ~1° scales are much less dramatic. Here, we will adopt an opti- 
mistic 1% systematic background uncertainty for dwarf galaxy observations. 

Relevant instrumental systematics are the spatial and spectral variations 
in the effective area, incomplete rejection of cosmic rays, and uncertainties 
in the energy and directional reconstruction (in case of the LAT, spectral 
uncertainties in the effect ive area range e.g. from 2% to 10%, depending on 



signature of interest |200| ). We will adopt here an optimistic reference value 
of 1% as uncertainty for the effective area, which still leaves room for future 
improvements in the instrumental design. In case of lACTs, we furthermore 
adopt the aggressive scenario that all hadronic showers are rejected, and 
only the cosmic-ray electron flux remains as an irreducible background (at 
energies below a few hundred GeV, this approximation is actually already 
realized with current technology). 

We will here consider two conceptionally different targets, which are of 
particular interest for current DM searches: a continuum signal from dwarf 
spheroidal galaxies and a gamma-ray line signal from the GC. 

Dwarf Spheroidals. We will here consider a single prototypic dwarf sphe- 
roidal galaxy; a combined analysis of several dwarfs could further improve 
the results if background and instrumental systematic are under control. We 
adopt a reference J-value of J = J dil J^^di p"^ = 10^^ GeV^ cm~^. inside 
an integration cone with radius 6 = 0.15° . This is of the order of what is 
expected for e.g. Draco or Segue I jssi lllol |: for other values our limits would 



roughly scale like oc J 9 . In case of space-based instruments, we estimate 



the level of the diffuse background by the IGRB determined by Fermi |147 . 
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For ground-bas ed in struments, we add on top of that the e^-flux measured 



by Fermi-LAT [201 



At the right end of Fig. [6|, with dark blue and red markers, we show the 
'fundamental' limits that result from the requirement that the signal flux 
from our reference dwarf is below a factor of 2a/2% (i.e. the 2a limit derived 
from two 1% errors) of the background within the integration region 9 at 
all energies. As apparent from this plot, especially space-based instruments 
still allow an impressive improvement of the limits if the effective area or 
observational time can be further increased. 

Gamma-ray lines. To estimate the fundamental limit for gamma-ray line 
searches in case of space-based instruments, we assume the same ROI and 



background fluxes as in Ref. [35|; the ROI has a size of ~ 20° and includes 
the GC. We adopt an energy resolution of 1% (as expected for e.g. GAMMA- 
400). Requiring that the peak of the line signal after convolution with the 
instrumental response (assuming a Gaussian for simplicity) does not over- 
shoot the background flux by more than 2a/2% yields the limits shown by 
the green markers at the right end of Fig. |6l These fundamental limits are 
10-100 times stronger than what is currently obtained with Fermi-LAT (see 
also Ref. jiof for prospects of observing line or IB features with current and 
future lACTs). 



7. Conclusions 

In this review, we have argued that one may consider gamma rays as the 
golden channel of indirect searches for DM in view of the extraordinarily rich 
spectral and angular information they can carry. This does not only help to 
discriminate signals from backgrounds but could eventually reveal valuable 
details about the properties of the DM particles. We have discussed the 
most important signatures in quite some detail and provided an update on 
current limits, demonstrating that indirect searches start to probe realistic 
cross sections and thus become competitive probes of physics beyond the 
standard model. 

While too early for a final judgement at the time of this writing, the line 
feature at 130 GeV that is seen in the Fermi data might turn out to be the 
most promising DM signal claimed so far. In fact, the intrinsic width of this 
feature must be smaller than roughly 20% - which leaves lines, VIB or box 
signals (Fig. [1]) as possible channels for an explanation in terms of DM. On 
the other hand, it is extremely challenging to find any explanation related 
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to astrophysics for such a spectral feature; for example, even a very hard 
contribution to the gamma-ray flux, with a sharp break at 130 GeV, cannot 
describe the data in a satisfactory way (Fig. [2]). The signal morphology is 
perfectly consistent with annihilating DM and an Einasto profile for the DM 
density, at least for distances larger than the apparent displacement from the 
GC by 1-2 degrees, but essentially rules out both cored and more contracted 
profiles (Fig. [3]); decaying DM is also in strong tension with the data. The 
data show a weak hint for a second peak at 114 GeV (Fig. [5] and Tab. [2]) 
which is exactly the combination of energies expected for the annihilation of 
DM particles into 77 and 7Z final states. However, large annihilation rates 
into these channels rather generically imply large annihilation rates rate also 
at tree-level, in potential conflict with continuum gamma-ray limits; VIB, on 
the other hand, does not suffer from this drawback (Fig. H]). 

If confirmed by the Fermi collaboration or other experiments, and in 
the absence of satisfactory instrumental or astrophysical explanations, this 
signal would lead to the exciting conclusion that the first particle beyond 
the standard model has been found in space rather than at a collider. We 
have discussed at length how astrophysical observations would already now 
help to determine detailed properties of this new particle. The situation 
will further improve in the relatively near future as prospects to study the 
130 GeV feature in more detail are extremely good. However, we believe that 
even if the DM origin of the signal is eventually not confirmed, our analysis 
serves to make a compelling case for the importance of focussing on clear 
spectral features in future searches for DM. 

During the last ten years or so, most limits on DM annihilation have im- 
proved by about one order of magnitude and this trend is expected to con- 
tinue for the next decade (Fig. [H]). We have also estimated the 'fundamental' 
(i.e. systematics-limited) reach of gamma-ray experiments, demonstrating 
that there is still quite some room for improvement even beyond those limits 
expected for the next decade, especially for space-based instruments. Even- 
tually, it may thus in principle be possible to probe cross sections down to at 
least one order of magnitude below the thermal value for TeV-scale particles; 
for many models, this would correspond to interactions too feeble to show up 
in any other kind of experiment, including direct or collider searches. While 
even those limits may not be sufficient to completely close the window for 
WIMP DM, model-building would need to become increasingly sophisticated 
to avoid them. 

Let us finally stress that in order to fully identify the properties of the 
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DM particles, it will of course be indispensable to correlate a suspected DM 
signal in gamma rays with results from indirect searches at other wavelengths 
and with other messengers. The same holds for direct searches and new data 
from colliders, both of which are guaranteed to deliver substantial new input 
in the near future — be it in terms of greatly improved limits or actual first 
hints for a signal. Chances are thus high that the next decade will either 
bring us a great deal closer to the long-sought nature of DM or, in the most 
pessimistic scenario in terms of detectional prospects, force us to seriously 
question the very idea of DM being composed of WIMPs. 
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